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, Abstract 

m ■ 

. Isoscaling and its relation to the symmetry energy in the fragmenta- 

\^ I tion of excited residues produced at relativistic energies were studied 

. in two experiments conducted at the GSI laboratory. The INDRA 

' multidetector has been used to detect and identify light particles and 

D , fragments with Z < 5 in collisions of ^^C on ^^^'^^^Sn at incident 

energies of 300 and 600 MeV per nucleon. Isoscaling is observed, and 
the deduced parameters decrease with increasing centrality. Symmetry 
term coefficients, deduced within the statistical description of isotopic 
scaling, are near 7 = 25 MeV for peripheral and 7 < 15 MeV for central 
collisions. 

, In a very recent experiment with the ALADIN spectrometer, the 

' possibility of using secondary beams for reaction studies at relativistic 

energies has been explored. Beams of lo^Sn, i24sn, i24La, and ^^Uu 
were used to investigate the mass and isospin dependence of projectile 
fragmentation at 600 MeV per nucleon. The decrease of the isoscal- 
ing parameters is confirmed and extended over the full fragmentation 
regime covered in these reactions. 
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1 Introduction 



The symmetry energy and its density dependence have received increasing 
attention in recent years because of their importance for nuclear structure 
and for astrophysics. Supernova simulations or neutron star models require 
inputs for the nuclear equation of state at extreme values of density and 
asymmetry 1J4 i2j . The strength of the symmetry term at about normal den- 
sity is quite well known from the masses of finite nuclei but its density depen- 
dence is experimentally very poorly constrained. Theoretical predictions are 
quite consistent for low-density nuclear matter .3] but diverge considerably 
at higher density where they are important for astrophysical phenomena as, 
e.g., the cooling of neutron stars |^. Nuclear reactions offer the possibility 
to investigate nuclear matter at densities other than normal nuclear den- 
sity, and an active search for observables suited to probe the strength of the 
symmetry term in reaction experiments is presently underway ^iSjiBj. 

Multifragmentation is generally considered a low-density phenomenon. 
However, the short-range nature of the nuclear forces causes a clustering in 
nuclear systems that have expanded as a result of an initial compression or 
heating in the course of a violent nuclear collision. Nuclear matter is not 
homogeneous under such conditions and the strength of the symmetry term 
is thus not easy to predict. In statistical multifragmentation models, it is 
assumed that normal-density fragments are distributed within an expanded 
volume. The density is only low on average, and standard liquid-drop pa- 
rameters are used to describe the nascent fragments including their isotopic 
degrees of freedom. In the Copenhagen version of this model (SMM), the 
symmetry energy term -Esym = — 2Z)'^/A is used with coefficients in the 
range 7 = 23 to 25 MeV [ZIIH]. 

The knowledge of the strength of the symmetry term under such con- 
ditions is of particular interest because the subnuclear densities and the 
temperatures reached at the freeze-out in fragmentation reactions overlap 
with those expected for the explosion stages of core-collapse supernovae [0]. 
This similarity permits laboratory studies of the properties of nuclear matter 
in the hot environment similar to the astrophysical situation. Nuclear mat- 
ter in equilibrium at low density is expected to consist of nuclei with a wide 
range of masses and with properties that dependent strongly on the global 
parameters of temperature, density, proton fraction and, in particular, also 
on the strength of the symmetry term under these conditions (Fig. ^ . 
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Figure 1: Mean charge-to-mass ratios (left top panel), and mass distributions 
of hot primary fragments (other panels) calculated with the SMM generalized for 
supernova conditions. Results are shown for T = 3 MeV and fixed lepton fraction 
Yl = 0.2 (left panels) and T = 1 MeV and fixed electron fractions = 0.4 and 
0.2 (right panels). Baryon densities are chosen as indicated. The dash-dotted lines 
represent calculations for density 0.1 po and a reduced symmetry term coefficient 
7 = 15 MeV (from 

2 Statistical description of isoscaling 

Isotopic scaling, also termed isoscaling, has been shown to be a phenomenon 
common to many different types of heavy ion reactions |S1 IIUI 1111 112j . It 
is observed by comparing product yields from otherwise identical reactions 
with isotopically different projectiles or targets, and it is constituted by an 
exponential dependence of the yield ratios R2i{N, Z) measured for the two 
reactions on the neutron number N and proton number Z of the produced 
fragments. The scaling expression 



R2i{N, Z) = Y2{N, Z)/Yi{N, Z) = C ■ exp{a ■ N + /3 ■ Z) (1) 



describes rather well the measured ratios for a wide range of complex parti- 
cles and light fragments 

In the grand-canonical approximation, assuming that the temperature 
T is about the same, the scaling parameters a and /? are proportional to 
the differences of the neutron and proton chemical potentials for the two 
systems, a = Afi^/T and (3 = A/ip/T. Of particular interest is the pro- 
portionality of the scaling parameters with the symmetry energy term. It 
has been obtained from the statistical interpretation of isoscaling within the 
SMM [S] and Expanding-Emitting-Source Model and confirmed by an 
analysis of reaction dynamics Jl]. The relation is 

A^,^ = ^A-A^ -47(|| -%) = -47A(ZV^^) (2) 

where Zi,Ai and Z2,A2 are the charges and mass numbers of the two sys- 
tems. The difference of the chemical potential depends essentially only on 
the coefficient 7 of the symmetry term and on the isotopic compositions. 

There is a simple physical explanation within the SMM why isoscaling 
should appear in finite systems. Charge distributions of fragments with fixed 
mass numbers A, as well as mass distributions for fixed Z, are approximately 
Gaussian with average values and variances which are connected with the 
temperature, the symmetry coefficient, and other parameters. The mean 
values depend on the total mass and charge of the systems, e.g. via the 
chemical potentials in the grand canonical approximation, while the vari- 
ances depend mainly on the physical conditions reached, the temperature, 
the density and possibly other variables. For example, the charge variance 
<7z ~ vT-^^/^t) obtained for fragments with a given mass number A in 
Ref. jl5j is only a function of the temperature and of the symmetry term 
coefficient since the Coulomb contribution is very small. 

The above relation opens the possibility for an experimental program for 
studying the evolution of the role of the symmetry term for the fragment 
formation as a function of the conditions of the experiment. Besides the 
isoscaling coefficient, the required inputs are the temperatures at breakup 
and the difference of the neutron-to-proton ratios N/Z of the two systems. 
An encouraging result was obtained when this method was first applied 
to the data of light-ion (p, d, a) induced reactions at relativistic energies 
of up to 15 GeV iSj. The obtained value had about standard magnitude 
7 23 MeV. This is not unexpected in this case because the data were 
inclusive. These reactions proceed mainly through heavy-residue formation, 
and the mean multiplicities of intermediate-mass fragments are correspond- 
ingly small (16j . 
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Multifragmentation becomes a dominant channel in reactions of heavier pro- 
jectiles or targets at relativistic energies. For the ^■^C on ^^^'^^^Sn reactions, 
studied with the INDRA multidetector pfS] in experiments performed at the 
GSI, maximum fragment production occurs at central impact parameters, 
according to the systematics ^Hl- The measurements were performed with 
enriched targets of "^Sn (98.9%) and ^^^Sn (99.9%). Reaction products 
with Z < 5 were detected and isotopically identified with high resolution 
by using the calibration telescopes of INDRA which are positioned at polar 
angles 45° < 6'iab < 176° 
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Figure 2: Scaled isotopic ratios S{N) for ^^C + "24243^ ^ 399 MeV (left 

panel) and 600 MeV (right panel) for three intervals of reduced impact parameter 
with "central" indicating 6/feniax < 0.4 and with offset factors of multiples of three. 
The dashed lines are the results of exponential fits according to Eq. Only 
statistical errors are displayed (from |20j). 

The ratios of the energy-integrated fragment yields measured for the two 
reactions and integrated over energy and angle obey the laws of isoscaling. 
This is illustrated in Fig. [21 for both incident energies and for three impact 
parameter bins, selected on the basis of the charged-particle multiplicity 
measured with the full detector. The fitted dependence on Z has been used 



to scale the ratios so as to obtain a single branch for each data set. The 
resulting slopes of the scaled ratios as a function of the neutron number N 
change considerably with impact parameter. The obtained fit parameters 
extend from a = 0.61 to values as low as a = 0.25 for the most central event 
group at 600 MeV per nucleon. 
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Figure 3: The left panels show the isoscaling coefficient a (top), the double- 
isotope temperatures TncLi (middle) and the resulting 7app (bottom) for the "'^^C on 
ii2,i24gj^ reactions at E/A = 300 MeV (fuU symbols) and 600 MeV (open symbols) 
as a function of the centrality parameter 1 — 6/6max- 

The top right panel shows the isoscaling coefficient a for hot (open circles) and cold 
fragments (dots) as a function of the symmetry term coefficient 7 as predicted by 
the Markov-chain calculations for ^^^'^^''Sn. The shaded area in the bottom right 
panel shows the region in the A(Z^/A^)-versus-7 plane that is consistent with the 
measured value a = 0.29 for central collisions and with the Markov-chain predic- 
tions for cold fragments. The dashed lines indicate the A(Z^/A^) = 0.0367 of 
"2,i24sn and the RBUU prediction (from EOI)- 

With temperature estimates obtained from the yields of ^'^He and ^'^Li 
isotopes, the values for an apparent symmetry term 7app, i.e. before se- 
quential decay corrections, were obtained according to Eq.[5J For peripheral 
collisions, they are close to the normal-density value 7 23 MeV but drop 
to lower values at the more central impact parameters (Fig. |21 bottom left). 
The effects of sequential decay for the symmetry term were calculated within 
the microcanonical Markov-chain version of the Statistical Multifragmenta- 



tion Model ^7] for excitation energies of 4, 6, and 8 MeV per nucleon and 
for 4 MeV < 7 < 25 MeV. The isoscaling coefficient a was determined from 
the calculated fragment yields before (hot fragments) and after (cold frag- 
ments) the sequential decay part of the code for which a standard value 7 
= 25 MeV was used. 

The hot fragments exhibit the linear relation of a with 7 as expected 
(Fig. 131 top right). For 7 smaller than 25 MeV, the sequential decay causes 
a narrowing of the initially broad isotope distributions which leads to an 
increase of the isoscaling coefficients a for cold fragments. The variation of 
a with 7 is thus considerably reduced with the effect that the value a < 0.3 
measured for the most central bins can only be reproduced with input values 
7 < 10 MeV. This is illustrated in Fig. 13 bottom right, which also shows the 
effect of possible variations of the isotopic composition of the two systems 
at breakup. Transport models predict that this difference should not change 
by more than a few percent j21j . 

4 Projectile fragmentation with ALADIN 

In a very recent experiment with the ALADIN spectrometer, the possibility 
of using secondary beams for reaction studies at relativistic energies has 
been explored 22 . Beams of ^■^^La, and ^^^Au were used to 

investigate the mass and isospin dependence of projectile fragmentation at 
600 MeV per nucleon. The neutron-poor radioactive projectiles ^''^Sn and 
^^^La were produced at the Fragment Separator FRS by fragmentation of a 
primary beam of ^^^Nd and delivered to the ALADIN experiment. Natural 
Sn targets with areal density 500 mg/cm^ were used in order to fully cover 
the rise and fall regimes of multifragmentation jl9j . A cross sectional view of 
the experimental setup and a short description are given in the contribution 
of De Napoli et al. to this workshop |23j . 

The mass resolution obtained for projectile fragments entering into the 
acceptance of the spectrometer is about 1.5% for fragments with Z > 6 and 
increases up to about 3% for the lightest fragments. The error in the mass 
determination has about equal contributions from the measurement of the 
magnetic rigidity {AR/R ~ 1%) and from the time-of-flight measurement 
(At ~ 150 ps, increasing up to 250 ps for light fragments). Masses are indi- 
vidually resolved for fragments with atomic number Z < 10. The elements 
are individually resolved over the full range of atomic numbers up to the 
projectile Z with a resolution of AZ < 0.2 obtained with the TP-MUSIC 
IV detector |22ll2SI. 
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Figure 4: Neutron multiplicities of the spectator source for ^''^'^^^Sn projectile 
fragmentation as a function of Abound (preliminary result). 

Results of an analysis of the Z fluctuations of the largest fragment within 
a partition are presented in f23; . The preliminary analysis of isotope ratios 
obtained for the studied reactions shows that isoscaling is observed. The 
decrease of the isoscaling parameters with increasing centrality is confirmed 
and extended over the larger fragmentation regime covered in the collisions 
with the Sn target. If normalized to the larger difference in N/Z that is 
achieved by using radioactive beams the new data are found to be consistent 
with the results obtained for ^■^C on ^^■^'^^^Sn. 

There is justified hope that the detection of projectile fragments without 
threshold with the ALADIN spectrometer will eventually permit the recon- 
struction of the spectator system at breakup. Fragments with Z > 3 exhibit 
well defined rapidity distributions centered near projectile rapidity 19 . It 
will be necessary, however, to also identify the spectator sources of hydrogen 
and helium ions and of neutrons. 

Neutrons emitted in directions close to ^lab = 0°, are detected with the 
Large-Area Neutron Detector (LAND) which covers about one half of the 
solid angle required for neutrons from the spectator decay. A preliminary 
analysis of the invariant multiplicity distributions of neutrons has led to the 
identification of the spectator sources of neutrons. They are characterized 
by temperatures up to about 4 MeV possibly caused by large contributions 
from evaporation. Reconstructed multiplicities for ^'^^■^^'^Sn projectile frag- 
mentation as a function of .Abound are shown in Fig. ^ On average about 3 
to 4 more neutrons are observed for the more neutron rich spectator system. 



5 Conclusion 

The present study of the role of the symmetry energy in spectator decay 
is motivated by the importance of this quantity for our understanding of 
astrophysical phenomena hke supernovae and neutron stars. The similar- 
ity of the thermodynamic conditions reached in nuclear multifragmentation 
and in supernova explosions gives experimental access to properties of hot 
supernova matter and its nuclear composition. The observed reduction of 
the symmetry term coefficient needed to describe isotopic distributions is 
consistent with a global classification of multifragmentation as a low-density 
phenomenon. More specifically, it may indicate that small and intermediate- 
mass nuclei are produced which are surrounded by other fragments and by a 
nucleon gas. The structure of these highly excited fragments, at the chemi- 
cal freeze-out point, may be considerably different from that of stable nuclei. 
The combined effects of the interaction with the neighbours and of their own 
expanded structure provide sufficient reason for a strong reduction of the 
symmetry term. 

In order to draw firm conclusions the evolution of the spectator sys- 
tems prior to breakup will have to be followed in more detail. The new 
experiments performed with the ALADIN spectrometer and with secondary 
beams offer the possibility to reconstruct the ncutron-to-proton ratios of 
the systems at breakup. A spectator source of neutrons has been identified 
by analyzing the coincident data measured with LAND. The secondary de- 
cay of excited fragments has the effect of reducing the widths of the isotope 
distributions and of narrowing their overall separation in N jZ . Model calcu- 
lations indicate that a reduction of the symmetry term will be partly masked 
by the sequential decay, with the effect that the true reduction at breakup 
may be considerably larger than what is observed for the final fragments. 
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